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Abstract 

In this paper we present a very simple and independent argument for the absence of the Boulware- 
Deser ghost in the recently proposed potentially ghost-free non-linear massive gravity. The limitation 
is that, in its simple form, the argument is, in a sense, non-constructive and less explicit than the 
standard approach. However, the formalism developed here may prove to be useful for discussing the 



, formal aspects of the theory. 

»hJ 1 Introduction 

It has been known for the very long time that giving a mass to the graviton in a stable and viable 
' manner is a very difficult task, if not impossible. At the linear level the only healthy mass term is the 

one found by Fierz and Pauli in their classical paper pQ. It propagates only five degrees of freedom in 
accordance with the general expectations for a massive spin-2 particle, while generically there would be 
six independent variables in the theory, with the sixth one representing a ghost. Unfortunately, the linear 
■ theory anyway contradicts observations due to the scalar graviton which couples to dust modifying the 

effective gravitational constant, but not to radiation keeping the bending of light intact. It was however 
later argued by Vainshtein [2] that non-linear effects will take over at small scales and save the whole 
day. But, almost at the same time, Boulware and Deser have shown [3] that even with the Fierz-Pauli 
mass term, be there any non- linear Vainshtein mechanism or not, the sixth degree of freedom comes back 
at the non-linear level reintroducing the ghost mode. And therefore a stable theory of massive gravity is 
probably not possible at all. 

However, recently a proposal for a ghost-free massive gravity has appeared [U [5J [5J. This theory 
has been extensively analysed in the perturbation theory, and the absence of the sixth mode was shown 
explicitly up to the fourth order in perturbations [71 [5]. At the same time, a fully non-linear Hamiltonian 
analysis has been done and proved the existence of the Hamiltonian constraint non-perturbatively [pj [TO] • 
After that, some doubts were expressed in the literature [IT] as to the existence of the secondary constraint 
generated by the Hamiltonian one which is needed for a consistent elimination of the sixth degree of 
freedom. But very recently the secondary constraint in non-linear massive gravity was evaluated (almost) 
explicitly [12] thus finalising the proof of absence of the Boulware-Deser ghost in this class of models. It 
justifies the large amont of interest which has been drawn towards understanding the phenomenological 
consequences of the new model. And now we have a rich literature discussing the black holes [TO1H^1[T51[TO] 
and cosmological solutions [T7J [TOJ [TO] in the massive gravity, as well as the first interesting results 
concerning the cosmological perturbations [2UJ . 

The ADM analysis of [9l [TOj [12] is technically quite involved; and given the paramount importance of 
the topic, we find it necessary to have a thorough understanding of this non-linear phenomenon of ghost 
exorcision from various vantage points. In this paper we would like to offer very simple arguments for the 
absence of the sixth degree of freedom in de Rham-Gabadadze-Tolley (dRGT) gravity at the fully non- 
perturbative level. In Section 2 we give a brief review of the standard ADM analysis and its application 
to massive gravity models. In Section 3 we introduce our set-up, and explain a very simple reason for 
the Hamiltonian constraint to exist, at least in the minimal dRGT-gravity with the fiat reference metric. 
In Section 4 we comment on the general dRGT models. And finally, in Section 5 we conclude. 



1 



2 Review of the ADM Hamiltonian analysis 



The standard way [2T] of doing the Hamiltonian analysis in GR is via the (3+l)-decomposition of space- 
time: 

ds 2 ee g^dx^dxv = — (N 2 — N k N k )dt 2 + 2N i dx i dt + lvj dx l dx j (1) 

where N and Ni are the lapse and shift functions respectively, and iVj ee jikN k . One can also find the 
inverse metric in terms of the lapse and shift functions and inverse spatial metric: 




ij _ TTNl 



(2) 



And as the Moo-minor of the matrix equals 7 ee detTy, we can conclude from g 00 = —jf? that 
\J—g = N-yff due to the standard rule of inverting the matrices. 

The next step is to calculate the Einstein-Hilbert Lagrangian density in terms of the ADM variables 
which is not so easy a task unless one takes the general geometric relations for embedded geometries 
directly from the textbooks. We find it most reasonable to recall the definition of the Riemann tensor 
as a commutator of covariant derivatives, and then to follow the geometric path outlined in the classical 
volume by Misner, Thorne and Wheeler [22]. But anyway, the result is 

y/—gR — Ny/j I R +K l k K k — (Kl) 2 I + (total time derivative and covariant divergence terms) 

where we have introduced the extrinsic curvatures 

1 /(») (3) \ 

Kik= 2N\ Vl Nk + Vfe Ni ~ ^ k J ' (3) 

and the three-dimensional scalar curvature and covariant derivatives as well as raising and lowering of 
the indices are defined by the spatial slice metric 7y . 

Now one can define the canonical momenta 7r u ee -J^- for the physical variables and also find the 
primary constraints ttn = and 7rjVj = for the momenta of the lapse and shift functions which act as 
Lagrange multipliers. The Hamiltonian then reads 



(3) 1/1/ ,\ 2 



H = - / d 3 x^ [ N ( R +- ( - (tt]) - n lk n lk ) ) + 2N l V " ) . 1 4 ) 



(3) 



The commutation of this Hamiltonian with the unphysical momenta directly gives the four physical 

((3) / 2 \ \ (3) 

R +i ( \ ~ n ik nlk J j an d Ci = —2^/7 y fc ^ifcj which reduce the number 

of propagating degrees of freedom from six to two, and make the Hamiltonian equal zero (in the weak 
sense of Dirac) as it should be in a time-reparametrization-invariant theory. 

2.1 dRGT gravity 

In the massive gravity theories one has to introduce an additional reference metric f^ v which can either be 
taken fixed (for example, Minkowski one) or endowed with its own dynamics thus producing a bigravity 
model, for otherwise there is no way to construct a non-derivative invariant which would act as a potential 
term for the graviton. (We will exclusively take the fixed f^ v metric option.) The basic ingredient of the 

dRGT model is the square-root matrix (\/g^f\ , and the minimal model potential is taken to be [B] 



V = 2m 2 



(v^T-a). (5) 



In what follows we will ignore the —6m 2 term which is introduced in order to avoid a contribution to the 
cosmological constant. And in the simplest case of Minkowski reference metric we deal with the square 
root of the following matrix: 



g^ aq (]av — N N j , - N N i N j ■ (6) 
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It is obvious now that the lapse and shift functions enter the Hamiltonian 



H = - J d 3 x^ (n (i? +i Q (^') 2 - 7r ife 7r ifc ^ -Vj+ 2N l m^j (7) 

non-linearly in the ^FfNV term, and therefore, naively one would expect to get some non-trivial equations 
for these unphysical variables instead of the physical constraints C and C l . If it was the case, we would end 
up with six degrees of freedom including the Boulware-Deser ghost. However, the very peculiar feature of 
the potential ((5} is that one combination of these constraints does survive as a physical relation leaving 
us with a healthy number of degrees of freedom [5] [8] . And we would like to understand the reasons 
for that. 

The standard approach to the Hamiltonian analysis [9j QUI [12] is to explicitly calculate the square 
root matrix in ([5]). Note that if it was not for the spatial metric 7, then the square root would have been 
really easy to find. Indeed, it follows from the very simple relation: 

1 a 1 \ ,„ fe k \ I 1 a 

-a? —a 1 a? 



( 1 " flV ) -at -aV ■ W 



However, one can not just simply take the square root of this part and then combine it with the square root 
of 7 _1 as they can never anticommute. But one can hope to redefine the shift functions Ni = (Sf +ND 3 i )nj 
such that the remaining A-dependent part of the \J g^ 1 / would still acquire this nice form [9] : 

V9~^f = — 7= , =f [ _„j _„i„j ) + ( n Yij(~ ) ' (9) 



Ny/1 - n k n k V ~ n3 ~ n * nJ J V xlJ {'J,n> 

and then the anticommutator with the X matrix would account for the difference between N l and n*, and 
between X 2 and 7. If this is so, then after such a redefinition the lapse function will enter the Hamiltonian 
linearly enforcing a truly physical constraint. Clearly, if there exists a combination of lapse and shifts 
which enters the Hamiltonian linearly then it must be possible to perform such a decomposition of \/g~ l f 
after a linear in N redefinition of shiftfl Indeed, in the (singular) limit of N — > 0, which is equivalent to 
rii —> Ni, the -^-part of the square root must tend to satisfying the property (|8|); and therefore it should 
always satisfy it because the difference between rij and Ni is determined by N on which no explicit 
dependence is allowed. And vice versa, if one finds such a change of variables, calculates the A 4 - 7 matrix 
and proves the relation ©, then the theory is free of the Boulware-Deser ghost. And it was actually 
done in [9] with the only limitation that the Dj operator is determined as a non-linear function of n, and 
therefore, in terms of initial variables, the transformation D is found only as an implicit function of the 
lapse and shifts. 

These results [5] have proven that the dRGT gravity is a potentially healthy deformation of GR at 
the full non-perturbative level. Moreover, it has been done for an arbitrary reference metric |10) and 
for non-minimal models too. The latter actually correspond to adding two more potential terms, V2 — 

(TrV^f-Tr (VF 1 /)' and V 3 - (Tr^T 1 /) 3 -3 (Tr^F 1 /) (tt (V^P/) 2 ) +2Tr (VfF^f '• 



3 The simple argument 

In our approach we introduce an extra matrix of auxiliary fields <§>£ into the model, so that the potential 
takes the following form: 

V = ! w(^ + ^~ 1 )> 2 9 ua u) (10) 

which yields the standard 2m 2 Tiy / g~ 1 f term <j5j) after integrating out the auxiliary fields. With the 
Minkowski reference metric we can safely demand $ k = $1 and <1>° = —& l , while in a general case 
some other combinations of $s will drop out of the actiono Now, if we set out to make every single 
step explicitly, the primary constraints are C\ = itn, Ch = t^n 1 an d — The constraint C3 will 

generate the matrix constraint C4 = <I> 2 — N 2 g _1 f ', and we are particularly interested in the constraints 
C5 and Cei generated by C\ and C2i respectively. 

1 If it was not linear then the lapse would definitely go non-linear in front of C's. 

2 If it does not seem so obvious as for how to arrive at the form mot , one may start with NV = 2m 2 «J>£I + 
K t ~ N2 9 va fav) and integrate out the KB. 
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With a simple commutation we obtain 



(S) 



where the derivative of the matrix ([6]) can be easily found to be 

T 2 <> «a, ( 5> 



dNl 9 lev y _ S k _ § J N k_ 5 k N j 

And therefore these constraints allow us to determine the shift functions in terms of 7 j&, ir lk and <£>; and 
a naive expectation would be that, in combination with C4, it is possible to express them solely in terms 
of 7^ and n lk . Then we find the last remaining constraint at this stage 

C 5 = + rn 2 V7 ^V%, = V7 (" S ~ Q (^) 2 - + 2m 2 iV 7 «) 

where one should not worry too much about summing over two upper indices as we just do not write out 
the unit matrix from the spatial part of the reference metric explicitly^ 
Up to that stage of analysis, the total Hamiltonian density is 

7 q\ 1 / 1 / \ 2 \ \ ( ,f;f ) / 



W = -VlN [ R +- I - (vrj ) - ir lk ir lk jj - 2^N< V k n lk + ^m 2 + (<S>-% N 2 g» a U) + 

+ VW ('-2v' ! 7r ifc +2m 2 (($- 1 ); + ($- 1 );iV^) > j , (11) 



and it would be its final form, with the set of purely second class constraints, for a generic choice 
of potential. But in the case at hand one can easily check that a particular combination of unphysical 
momenta (see the subsection below) does actually commute, in the weak sense, with the total Hamiltonian 
irrespective of the values of Lagrange multipliers. Hence, the constraints C4, C5 and C§ do not allow to 
unambiguously express the naive unphysical variables, N, Ni and $, in terms of the spatial metric and 
its momenta. And we can nothing but conclude that they do contain a non-trivial constraining equation 
for the would-be-physical variables 7^ and n 1 ^ Q This observation is actually the final step in our proof 
that the dRGT gravity contains strictly less than six degrees of freedom. 



3.1 Some technical details 

We now proceed to explicitly calculate the commutators of the unphysical momenta with the total 
Hamiltonian. Obviously, within the constraint surface, we only need to commute them with the other 
constraints because the commutations with the first line in the Hamiltonian (JTTJ) have already been done, 
and the very meaning of the other constraints is that those commutators do weakly vanish. Therefore, 
we can find 

-±={Ci , H} = -2A 4 J fe AY fc + 2A 5 m 2 t^ 1 )] 7 « 
-±-{C 2t , H} = -2A4 + 2X 4 lN k + 2m 2 K k , 



3 Note also that we could define $ 2 = g~ 1 f instead of <J> 2 = N 2 g —1 f, and as one can easily show, it would have resulted 
in the contribution to the constraint C5 which could be brought to a very simple form, 2m 2 ($ — 1 )* 7 lJ , by use of the 
C4 = equation. This form contains neither the lapse nor shifts, although the lapse would have appeared in C$i. And it 
would have been the first instance to face a crash of our naive expectations, that is, modulo some mixing, there are roughly 
three systems of second class constraints: C± and C5, C2 and Cq, and C3 and C4. However, despite its apparently striking 
form, this fact is not that easy in being promoted to an actual proof, and therefore we will proceed with the definition l|10jl 
which allows for more pleasant calculations. 

4 This is paralleled by the fact that in the standard approach the total Hamiltonian commutes, at this stage of analysis, 
with 7rjv with no restrictions on the values of Lagrange multipliers, and indeed we have a physical constraint instead of an 
equation for N . 
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-L{C 3 » , H} = 2\^Z + 2m*X 5 ^N^ (^)\ + 2m 2 V ^ + rf) 

where we use the symmetry of the auxiliary fields and their Lagrange multipliers, = A4JJ, and A^ = 
— A4q- Normally, these commutators would give us fourteen independent linear equations for fourteen 
Lagrange multipliers, so that all of them would be set to zero. However, using the C4 constraint and a 
simple formula ( < i 5 ~ 1 )^ -^s = — (3> -2 )^ (^ _1 )k> one can deduce from the last commutator the 

following relations: 

^ - ^= (*- x )l W ,H ]+m * + ( ( *-i)< + , 

And we see that the combination of A4J — X^N^ depends only on the commutators with C3, while both 
A5 and Xq completely drop out of this expression. And therefore, using the commutator with C2, we 
obtain 

2m 2 A 6 fc (<S>~% = ± {c 2l + ~ (* , h} 

in the weak sense. (Note that one should not worry about commuting the Hamiltonian with the coefficients 
in front of the constraints because the momenta do vanish on the constraint surface.) On the other hand, 
one can compute XijY J an d compare it with the {C\ , H} commutator. It (weakly) determines a linear 
combination of Ags: 

2m 2 A 6 ' £ ((<i>- 1 ) fc + ($- 1 ) i fc iV i ) =-^{N 2 (Z-y^Csf+NCr , H } . 
These two results must agree, and it singles out a combination of momenta 

TT N +N($ 1 ) a7 ^+A _ I I* J (^ + ( $ )„,**%-(* ) a N3 *9i) 

which weakly commutes with the Hamiltonian for any values of the Lagrange multipliers. 

In our approach, this combination determines the direction in the space of unphysical variables along 
which there has been no restriction so far, under any of the constraints. This corresponds to the indepen- 
dence of N in the standard treatment. In either approach, a constraint which commutes with all the other 
constraints (and leaves one combination of the Lagrange multipliers undetermined) may appear in two 
distinct situations: either it is a geniun first class constraint and corresponds to a gauge freedom in the 
model, or some extra constraints are needed for the self-consistency so that the whole set of constraints 
is non-degenerate second class. As the former seems not to be the case, the generation of a one more 
constraint is unavoidable (and, at the end of the day, the values of all non-dynamical fields should be 
somehow determined unless there is a gauge freedom indeed), and therefore the scenario forseen in the 
reference is a priori highly implausible. 

Technically, what follows at the next step is that the Lagrange multipliers Ai, A2, A3 should serve to 
ensure the preservation of C4, C5, Cq constraints. But a factor in front of the all-commuting combination 
of momenta drops out of this game; and there would not be enough freedom to make all the necessary 
commutators weakly vanish, if there is no gauge freedom. (And in our case the commutator with the 
C5 constraint has two independent parts, proportional to the delta-function and to its derivative, and 
the one undetermined combination of the other Lagrange multipliers would be used for fixing the latter 
part.) This is how an extra constraint is generated in the model. However, it is not a priori clear that 
it would be a physical one making the number of degrees of freedom not more than the healthy amount 
of five (as opposed to five and a half) because the next constraint could just add the fourth equation 
for the lapse and shift functions without any information on the spatial sector. Due to this reason it is 
harder to proof that the theory contains five degrees of freedom than just to show that its number is 
less than six. But a plausible argument would be that, after integrating out all but one undetermined 
unphysical variable, the remaining combination should go linearly in the action with the coefficient equal 
to the remaining constraint which would therefore commute with this part of the Hamiltonian (this is 
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subject to the criticism in [H]). And as is known by now [12], the strange 5 ^-situation is not the case, 
and two extra constraints do appear in the model, one is the "secondary" one for the spatial variables, 
and the other finally fixes the lapse and shifts. 

We would not proceed with explicit derivations in this paper because the most important result is 
already known [9j [12] , and our only purpose was to present an alternative and fairly simple method of 
analysis. However, all the necessary calculations are very straightforward although time consuming. 



4 On arbitrary reference metrics 
and non-minimal dRGT models 

Fow the sake of simplicity, up to now we have considered only the simplest choice of the reference metric, 
i.e. the Minkowski one. However, the non-linear massive gravity has been proven to be free of the 
Boulware-Deser ghost for any choice of the reference metric |10j . and even in its bigravity version too 
[T2] . Incorporation of an arbitrary lapse and an arbitrary spatial metric is actually trivial (and the latter 
even makes the location of the spatial indices nicer) , while shifts do produce some problems and affect 
the simple form of the decomposition ([9]) , see [10] . It can be readily seen by taking the general reference 
metric 

^(-(M^M*) (12) 
and calculating the basic building block of the model: 

M 2 -M k (M k -N k ) .s,,.V' A/, 



_V LL + 7 »J Mj s ik -f 



However, the complication is a relatively mild one: the C5 constraint receives a more involved con- 
tribution of 2m 2 N y /j ((r 1 )"]^ + (* _1 ) I fe Sfci) 7 ij instead of the simple 2m 2 N y /j l ij , and the 
form of Cq is also changed in an obvious way. Nevertheless, it simply corresponds to a rotation of the 
variables, and the subsequent calculations become a bit more complicated only due to the form of the 
coefficients with no crucial change to the results. It actually should have been the case because, at least 
in the class of coordinate-independent reference metrics, a general metric can be transformed to zero 
shifts by a (linear) change of coordinates, e.g by the one which diagonalizes the matrix /0 

The models with general potentially ghost-free potentials can, in principle, be treated in the same 
wa yll Indeed, in order to check that the constraints Ci, . . . ,Cq are preserved during the evolution, one 
has to perform a straightforward computation of the commutators and to decide upon the solvability of a 
system of linear equations for Ai, . . ., Ag- And if some more constraints are required then they definitely 
convey a non-trivial piece of information about the spatial sector of the model. 

However, once we have understood the reasons for the sixth mode to be absent in the minimal dRGT 
model, the ghost-free nature of the higher potentials can be better explained with the standard argument 
of the references [9j [TO]. We know that, after a linear in N change of variables, the minimal model 
action contains the lapse function only linearly. And if one has convinced himself that it implies the 
decomposition ([9J, then it is also obvious [9j[T0J that the same is true of the V2 and V3 potentials because 
for a matrix A with the property © we have Tr(^4") = (Tr^)", and all unwanted powers of the lapse in 
the potential do cancel. This is how the symmetric polynomials of the eigenvalues come into play; and 
in four dimensions there is only one more of them, det(g -1 /), which being multiplied by \J—g produces 
nothing but a constant shift of the action. Our method does not respect the decomposition ([5]), and 
therefore it is not so elegant in generalizing to non-minimal models. 



5 Note that one can, of course, make a coordinate transformation, even in massive gravity, as long as both the physical 
and the reference metric are being changed accordingly. And recall that the very property of *j g~ 1 f to be decomposable 
into the sum of N' 1 and N° parts can be proven at each point of the space-time manifold separately, with no reference to 
its coordinate dependence or independence. 

6 There is a somewhat subtle case of the purely quadratic potential V2 for which the variation of the <J>^ ($ — 1 )^ (9 1 /) " — 
(g -1 /)^ term with respect to $ would leave an overall scalar factor undetermined much like the conformal invariance appears 
in the Polyakov action of the bosonic string. 
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5 Conclusions 



We have presented a new method of non-perturbative analysis of the non-linear massive gravity. And in 
particular, we give a very simple argument for the absence of the Boulware-Deser ghost in this theory. 
The limitation of our approach is that the argument is non-constructive, and does not proceed directly 
in terms of the metric components. However, in principle, one can make all the derivations explicitly, 
and calculate the number of constraints and independent degrees of freedom. In general, the power of 
this approach is in the fact that many things can be done at the level of a bit bulky but absolutely 
straightforward calculations, very automatically, with no need of making more qualified and creative jobs 
such as taking the square root of a matrix. 

Admittedly, the introduction of ten more configuration space dimensions (and ten more pairs of 
constraints to eliminate them) is not very helpful for the actual physical calculations. But, given the 
relative ease with which this argument shows the presence of extra constraints in the model, it is reasonable 
to hope that it can be useful for discussing the formal aspects of the theory. And anyway, for the real 
calculations a somewhat different language is better suited [51 [71 [5]. Obviously, it would be very interesting 
to find out whether any kind of such formal tricks with auxiliary fields like $ could produce a similarly 
simple argument for the full non-linear stability in order-by-order perturbation theory. 
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